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Spin transport properties of a coupled bilayer electron gas with Rashba spin-orbit coupling are 
studied. The definition of the spin currents in each layer as well as the corresponding continuity- 
like equations in the bilayer system are given. The curves of the spin Hall conductivities obtained 
in each layer exhibit sharp cusps around a particular value of the tunnelling strength and the 
conductivities undergo sign changes across this point. Our investigation on the impurity effect 
manifests that an arbitrarily small concentration of nonmagnetic impurities does not suppress the 
spin Hall conductivity to zero in the bilayer system. Based on these features, an experimental 
scheme is suggested to detect a magnification of the spin Hall effect. 
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I. INTRODUCTION 

Manipulating the spin degree of freedom for electrons 
has recently brought in an emerging information technol- 
ogy, spintronics 0, H, S] , which offers novel clues for de- 
signing devices based on traditional materials with spin- 
related effects. In this promising field, the spin Hall ef- 
fect [3j S @] is regarded as a candidate method to inject 
spin current in semiconductors. Based on the spin-orbit 
coupling (SOC), an external electric field is required to 
drive a transverse spin current while the magnetic field 
is not necessary, which is much different from the tradi- 
tional applications of the spin degree of freedom. A uni- 
versal spin Hall conductivity c/Stt is predicted theoreti- 
cally in a clean single layer electron system Q. Several 
groups' calculations 0, S, 0, [13, Ell showed that non- 
magnetic impurities would suppress this spin Hall con- 
ductivity to zero while others indicated that the spin Hall 
conductivity is not zero in the presence of magnetic impu- 
rities [isj . Experimentally, the spin accumulation in 
nonmagnetic semiconductors has been observed (l^ . [isj 
and the spin current was detected either by Kerr rotation 
microscopy llal or by two-color optical coherence control 
techniques [17| . Very recently, a direct electronic mea- 
surement of the spin Hall effect has been reported [l^l 
where the spin current induces the charge imbalance and 
a voltage is detected. 

As the SOC, which is crucial to the spin Hall effect, is 
a relativistic effect and thus comparably weak, a natu- 
ral question is how to strengthen this effect. In the light 
of single layer systems being considered in current lit- 
erature, one may ask whether a multi-layer system pos- 
sesses a magnification effect and what new phenomena 
will take place if the tunnelling between layers is taken 
into account. Another more realistic question is what will 
happen if there exist impurities in a multi-layer system. 

In this paper, we investigate the spin transport prop- 
erties in a coupled bilayer electron system with different 
SOC strengthes in each layer as well as the tunnelling be- 
tween layers. As a starting point, we generalize the defini- 
tions of spin currents to a coupled bilayer system and ob- 



tain the corresponding "continuity-like" equations. Car- 
rying out calculations of the spin current in the Heisen- 
berg representation, we find that the spin Hall conductiv- 
ity in each layer manifests abrupt enhancement around a 
particular value of the tunnelling strength between lay- 
ers and undergoes a sign change across this point. The 
influence of impurities is also studied. We indicate that 
the spin Hall conductivity in the bilayer system can not 
be suppressed to zero by an arbitrarily small concentra- 
tion of impurities. An experimental scheme is designed 
on the basis of these features to magnify the spin Hall 
effect near the turning point. Besides, possible logical 
gates are expected to be elaborated based on the sign 
change of the spin current across this point. 

The whole paper is organized as follows. In Sec. [Hi we 
generalize the definition of the spin current in each layer 
and obtain the "continuity-like" equations. In Sec. IIIIl 
the spin current as well as the spin Hall conductivity in 
each layer are calculated in Heisenberg representation. In 
Sec. IIVI the influence of disorderly distributed nonmag- 
netic impurities on the spin Hall conductivity is investi- 
gated. In Sec.|Vl we show the greatly enhanced spin cur- 
rents near the turning point and scheme out possible ex- 
periment to detect a magnification of the spin Hall effect. 
Finally, a brief summary is given in Sec. IVII and some 
concrete expressions are written out in the appendix. 



II. "CONTINUITY-LIKE" EQUATIONS 

As a proposition to study the spin transport, we firstly 
introduce the definition of the spin current in a coupled 
bilayer system in this section. Throughout the whole 
paper, we consider a coupled bilayer system where the 
strengthes of the Rashba-type SOC in each layer are 
different and the tunnelling between layers always oc- 
curs. The spaces spanning the electrons' spin states and 
layer occupations, respectively, carry out SU(2) represen- 
tations. If the spin and layer representations are denoted 
by Pauli matrices a a and r-matrices r^, respectively, the 
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total Hamiltoiiian of such a system can be written as 



presence of SU(2) gauge potentials: 
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where ai and a2 refer to SOC strengthes in the front 
and back layers, correspondingly, and /3 the tunnelling 
strength between layers. / stands for the unit matrix. 
For convenience, a+ = (ai+a2)/2 and a_ = (ai — a2)/2 
are introduced in the second line of the above equa- 
tion. Hereafter, indices a and i run from 1 to 3. Let 
V'f = (</'f|j 0fx)"^ and f/'jj = (</*bt' ^bj)"^ represent the spin 
states of the electrons in the front and back layers, respec- 
tively. Hereafter, the layer-index f or b labels either the 
front or back layer. Then a four-component wave func- 
tion, denoted by = (0fp ^,3^, 0^,^)^ = (V'f,'0b)^ 
must be introduced for a complete quantum mechanical 
description of the system. 

The well accepted definitions of the spin density and 
the spin current density in a single-layer system are 
S"" = -^^s"-^ and J" = RevI'tj'Jvl/^ respectively. Here 
s° — (Jah/2 is the spin operator and = i{v, s"} the 
spin current operator with the curl bracket denoting the 
anti-commutator and v = 7^ [f > ^^o] the velocity opera- 
tor. The bold face manifests the quantity is a vector in 
the spatial space, e.g. J"^ — (J^, J" , J^). It is natural to 
define the full spin current operator for the whole bilayer 
system as 
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with j° and j{J being the spin current operators in the 
corresponding layers. Even though the tunnelling couples 
two layers, the spin current operator is in a block diago- 
nal form since the tunnelling is momentum- independent. 
Then we have the spin density and the spin current den- 
sity in each layer 



(3) 



where I stands for f or b. 

It is obvious that the presence of the SOC, which 
can be regarded as certain SU(2) gauge potentials Ai 
and Aq [l9l |. leads to the non-conservation of the spin 
density. Hereafter, a vector in the spin space is de- 
noted by an overhead arrow, e.g. S = {3^,3^,3^). In 
terms of these gauge potentials, the partially conserved 
spin current takes a covariant form and obeys the 

"continuity- like" equation, namely, (-^ — rjAox)3 + 



ox. 



0. Through an analogous procedure 



as in Ref. [19|, we can derive a general "continuity-like" 
equation for the spin density in each single layer in the 
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(i^ls^^-iPis^Pi). (4) 



In the coupled bilayer electron gas with Rashba SOC, 
Ax = |?(0, ai ,0), Aiy = ,0), Ax = 

fi(0, a2 ,0), Ay = -^(«2, ,0) and A. = Ao = 

Ahz—Aha = ^ with rj = h. The tunnelling between layers 
gives rise to the term on the right hand side of Eqs. (|4]) 
and this term results in additional non-conservations for 
the spin density in each layer. 



III. SPIN CURRENTS IN A CLEAN SYSTEM 

In this section, we calculate the spin currents for a 
clear system in Heisenberg representation po| . A weak 
electric field E = Ex applied on both layers is regarded 
as a perturbation. We mainly focus on J^^ component of 
the spin current in the €-layer which is flowing perpen- 
dicularly to the electric field with the spin polarized in 
the z-direction. 

Diagonalizing the unperturbed Hamiltonian ([1]), we 
obtain four energy bands: 
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/92 -f- a^k? — a+k, 



(5) 



sgn(a;) = 



1 if x > 0, 
if x = 0, 
-1 if a; < 0, 



and ki 



(3/ — o?_ denoting a special point where 

El = £2 = f?k\l2m. The landscape of these bands are 
plotted in Fig. [1] in which "f in the right panel marks the 
level crossing point of ei and £2 at kt. As we will see later, 
the spin Hall conductivity exhibits sharp cusps around 
this point. In the following, we consider the case k < kt 
which has the same result as fc > kt. The eigenvectors 
= (V'fj: V'bj)"^ with j — 1,2,3,4 labelling the band 
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FIG. 1: (color online) The four energy bands corresponding to 
the four eigenstates given in Eq. (|SJ . The surface of revolution 
in the left panel is obtained by revolving the curves in the 
right panel with respect to the vertical axis. The | in the 
right panel marks the level crossing point of ei and £2. 



with kox and koy being the initial values and 



— {l®(Tz) = - [a+k^I (g) ax + a+kyi (g) ay 

+a^kxTz® a^ + a^kyT^® ay\. (7) 



Obviously, the time evolution oi I ® az depends on those 
of other four-by-four Hermitian matrices, such as, I ® ax 
which also depends on other matrices. Hence, we need 
to deal with the time evolutions of sixteen matrices {I ® 

I,I(S)ax,I(E)ay,I(E)az,Tx(S)I,Tx(E)ax,Tx(E)ay,Tx(S)az,Ty(E) 

I, Ty Ig) ax,Ty ay, Ty ig) a z , T ^ I, T ^ g a X , T ^ g a y , T ^ g 

az} which span the space of the four- by-four Hermitian 
matrices. If we arrange those 16 matrices successively in 
a single column, denoted by F, the problem reduces to 
search solutions of a set of 16 linear differential equations: 



indices are given by 
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where tp = tan ^(ky/kx) and the normalization coeffi- 
cients Nj are given in the appendix. 

The spin current operator for the whole bilayer system 
1 h^k 
is given by jy = ^i^y ^ ^ ^"^} ~ 2 ^ ^ ® '^^^ 

lutions of operators are governed by Heisenberg equation 

eEt 

of motion. Thus we have kx = kox ^ and ky = koy 



2 / eEt 

5.r.-(M + -M.)r, 
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where the concrete expressions of M and Mt are given in 
the appendix. 

Expanding F in scries of the electric field, namely, F — 
p(o) _|_ p(i) -|- . . . ^ we have the following equations: 



5iF(°) = |mf("), 
n 

a^F^i) = |AfF(i) + ^MtF(") , (9) 



up to the first order. Using the standard method to solve 
these equations, we obtain the linear order term / (g) ai^' 
in the limit f — s- 0: 



eE 



(^Ci/(g) aox + C2l'S) aQy + C3Tz<S) aox + CiTz® croy) , (10) 



where aoa stand for the initial values of ctq at t = and 
the coefficients C are written out in the appendix. 

The spin currents in both layers produced by the states 
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in each energy band are evaluated as 



■) = -eEh^ sinV x 



+ a^fc2 + a_fc)[/32 - (a2_ -a+a_)fc2] 
8TOa+ fc y'/32 + a?.fc2(/32 - (a2_ - ) fc2 ) 
(V-i.bliyl^Ai.b) = -eEh^sin^ip x 

i^P^ + alk^ + a-k)[{al + a+a^)^ - /?2] 
8TOa+ fc y^/32 + a^fc2(/32 _ (a2_ - a?. ) fc2 ) 

(V'3,bby|V'3,b) = e^'^^ sin295 X 

(y^/?2 + a2_fc2-a_fc)[(a2_ + a+a_)fc2 - /32] 



while 



(V'2,t(b)|j^|V'2,t(b)> 
(V'4,f(b)|j^|V'4,f(b)) 



-('/'l,t(b)|j^|V'l,f(b)>> 
= -(V'3.f(b)|jy|l/'3,f(b)>- 



(11) 



The total spin current in each layer is the sum of the 
contributions of the four bands up to the Fermi level, 
J!(i,)y = Ej,fc('/'i,f(b)|jylV'j,f(b)>"Ffe)/(i:«; X Ly) 
where np{ej) is the Fermi distribution function and 
Lx X Ly the size of the system. The explicit expressions 
for "/f^^b) y at zero temperature are given in the appendix. 
Eqs. ([TT|) tell us that the spin currents produced by the 
states in bands £1 and £2 are always with the opposite 
sign. Thus only the contributions by the states in 62 with 
momentum kpi < k < kp2 remain. The case for bands 
£3 and £4 is similar. Here and throughout the paper, kpj 
denotes the Fermi wave vector in the band Sj. 

The full spin current of the whole bilayer system is 
given by Jy = Jfy + Jhy the corresponding spin Hall 
conductivity is defined as as = dJy/dE. Our results 
can also be verified by Kubo formula. It is worthwhile 
to observe our results in two specific cases. In the case 
that the tunnelling is absent, the system becomes a de- 
coupled two single layers and its spin Hall conductivity 
becomes as = e/47r, twice of the universal value in a 
single layer Q. In the case of a_ 0, there is no dif- 
ference between the two layers and thus they can not be 
distinguished. Consequently, no matter the tunnelling is 
present or not, they behave just like decoupled two sin- 
gle layers since tunnelling to the other layer makes no 
difference from staying in the original one. 

Now we are in the position to investigate the tunnelling 
dependence of spin Hall conductivities erf and Ub in each 
layer. Based on the above results, we plot erf, (Tb as well 
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FIG. 2: (color online) (a): Spin currents in each layer 
and in the whole system versus a_ are plotted with a+ = 
0.55 X 10"^^eVm and /? = 4.013 x 10"*eV. (b) and (c): Spin 
currents in each layer versus the tunnelling strength /3 are 
plotted with ai = 10~^''eVm and 02 = 10~^''eVm, in which 
sharp peaks emerge near Pt . The choices of the other param- 
eters in the above figures: the Fermi energy ep = O.leV and 
the effective mass m = 0.05me. 



as as in Fig. [21 The dependence of the spin Hall conduc- 
tivity in each layer on the strength of the SOC is quite 
different from that in a single-layer system which does not 
vary as the strength of the SOC changes. As illustrated 
in Fig. ma), a^^^-^ increases (decreases) monotonously as 
the strength of the SOC in its layer increases (decreases) 
while CTs keeps constant. We also plot <7f and a^ ver- 
sus the tunnelling strength in Fig. Hfb-c). They change 
abruptly near j3t where kpi = kp2 = kt and also undergo 
sign changes across this point. At this point, the spin 
current produced by the states in band £1 and that in 
band £2 cancels each other precisely, leading to a depres- 
sion of as which always keeps a constant value e/47r for 
(3 ^ (it- It manifests that each layer posses a large spin 
conductivity near Pt while as of the whole system re- 
mains e/47r. These features are instructive for designing 
experiments to detect a magnified spin Hall effect. 



IV. THE INFLUENCE OF IMPURITIES 

In realistic systems, disorderly distributed impurities 
are unavoidable, which frequently affect transport prop- 
erties. It was believed that the spin Hall conductivity 



in a single layer electron gas could be suppressed to zero 
by the vertex corrections of nonmagnetic impurities even 
for an infinitesimally small concentration [7l.l8l.[9l.[l0l.[lll|. 
Therefore it is necessary to investigate the impurity ef- 
fects in the bilayer system. 

We consider disorderly deployed nonmagnetic impuri- 
ties. The short-ranged interaction between the electron 
and impurities at positions Ri is described by Vim = 
^jU5(r — Ri). Here we assume the coupling strength 
u is sufficiently weak so that the Born approximation is 
applicable. The averaged retarded Green's function sat- 
isfies the Dyson equation = Gq + G^T^^G^ where 
the overline refers to an average taken over the configura- 
tion of impurities, G^ denotes the free Green's function 
and the self energy brought about by the impurities. 
In the Born approximation the Dyson equation can 
be explicitly written as 



Y^G^{q,u))G''{p,u)m 



G^{p,Lo)^G^{p,w) 



where n^j^ stands for the impurity concentration and v 
the size of the whole system. 

For convenience, we introduce the chiral representation 
in which the Hq is diagonalized by a unitary matrix U, 
i.e. U^HqU = diag(ei, Esj £4)- In this representation, 
the free retarded Green's function reads G^^^i_^^{p, uj) — 

diag(^(tj-£i+«r])"\ (tj-£2 + «7?)"\ ('^-£3 + «?7)"\ (w- 

£4 + irj)^^^ such that equation (fT2|) solves 
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with gj = l/{u; — Ej + ^) for j = 1, • • • , 4. Here t = 
{2TTu'^n-^^Np)^^ is the momentum-relaxation time and 
Np the density of states of the electron at the Fermi 
surface. 

In terms of the Kubo formula, the averaged spin Hall 
conductivity at zero temperature can be calculated. 



as (uj) 



-UJ 



X [ G^iuj'+cu) - G^{uj'+uj) ] je G^(uj') 

+ 0{-Uj') G^iLU'+L0)U G^iLO')~G^LO') ]},(13) 

where G"^ is the advanced Green's function, = ev^ 
the charge-current operator and 0{lu) the step function 
representing the Fermi distribution function at zero tem- 
perature. The trace Tr in Eq. (|13p implies both the con- 
ventional trace over the spin indices and the summation 
over the momenta. In the uncrossing approximation , 
as is the sum of ct^ and erf" , the former is the contribution 
by one-loop diagram 




while the later is that by a series of ladder diagrams 

qA qA 




je + j' 




Je + 



A. One-loop diagram contribution 

To derive the dc conductivity, we take the limit w — *■ 
in Eq. p3p and obtain the one-loop diagram contribution 

= :^Tr[j^ip)G^ip)Up)G^ip)] 
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.fci 



where x = a a" ^ function of /?/(«- 
£i{kp) — ej{kp) is the energy splitting between two bands 
at the Fermi surface. The Fermi wave vector fci? is given 
by ^i/2mfi/h with fi being the chemical potential. In car- 
rying out the summation of momentum in Eq. (|14p . we 
have adopted the large Fermi-circle limit 1/r, Aij . 

Our result in Eq. seems to be similar to the ex- 
pression of a single layer system, a^f^ = ~ i-t-A^r^ )• 
Moreover, the extra term in the last line of Eq. (fT4l) 



and the pre-factors ^re peculiar in the bilayer sys- 
tem. It is worthwhile to observe the aforementioned two 
specific cases. In the case of zero tunnelling /3 ^ 
ix = -1), we have (7° = ^(1 - + 1 - j^^) 

where A/ = 2aikF and A;, = 2a2kF are the spin-orbit 
splittings in each layer. It demonstrates that the sys- 
tem reduces to a decoupled one. In the twin-layer case 
Q!_ ^ (x = 0), we have a° = £-{1 - j^r^rfj) which is 
just twice of the value of a single layer system. This is 
actually a trivial case as there is no difference between 
layers even though the tunnelling is present. The above 
reasonable conclusions are consistent with the results for 
the clean system derived in previous section. 



B. Vertex correction 

The sum of ladder diagrams gives rise to a^. By in- 
troducing a matrix-valued vertex Jy which is the sum of 
the vertex correction to jy, diagrammatically, 
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then (j^ can be written as As a result, we have 

= ^T^A Jy G^iP) Up) G^iP) I (15) 

where Jy is momentum-independent and satisfies the 
transfer matrix equation 

Jy = ^T.^^('i)0'y{q) + J^)G''iq). (16) 

I 
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where J12, J14, J34 are the matrix elements of Jy. Their 
explicit expressions, a solution of Eq. p6|) . are given in 
Eq. (|A7p in the appendix. When the tunnelling vanishes, 



we have at' — — (1 ^ 



1 



l+Air- 



■), reducing to 



the case of a decoupled bilayer system, and precisely 
cancels , leading to a vanishing spin Hall conductivity. 
The nontrivial situation is that both the tunnelling (3 and 
the difference in Rashba strengthes a_ are present, which 
makes ct, survives. 



V. MAGNIFICATION EFFECT AND POSSIBLE 
EXPERIMENTS 

The spin Hall conductivity is the sum of ct^ and . 
An arbitrarily small concentration of nonmagnetic im- 
purities can not suppress the spin Hall conductivity in 
a bilayer system to zero, which is quite different from 
the case in the single layer system. In Fig. we plot 
the spin Hall conductivities for each layer (Tf(b) and for 
the whole system as with different parameters. Panel 
(a) of Fig. © is the plot for a+ = 0.55 x IQ-^^Ym 
and a- — 0.45 x 10~^^eVm {i.e. The strengthes of the 
Rashba spin-orbit coupling in each layer are of the same 
order). The curves for the conductivities exhibit sim- 
ilar cusps around the turning point Pt as in Fig. ([2|) 
without impurities. The conductivity in each layer pos- 
sess opposite signs, leading to a quite small as for the 
whole system. However, things are changed when a_ is 
comparably large. Panel (b) in Fig. ([3|) shows the con- 
ductivities with parameters a+ = 0.55 x 10^^'^eVm and 
a_ = 0.45 X 10~^'^eVm, i.e. the Rashba strength in the 
front layer is ten times as much as that in the back layer. 
The opposite signs of the conductivities in each layer in 
the absence of impurities turn to be the same in presence 



of impurities. As a result, the peak value of as for the 
whole system is considerably large. It suggests that a 
large difference in the strength of Rashba spin-orbit cou- 
pling between layers is favorable for a greatly enhanced 
spin Hall conductivity. 

Above results are obtained with rather dilute impuri- 
ties which requires the mobility of the two-dimensional 
electron gas to be quite high. The influence of the con- 
centration of impurities on the spin Hall conductivities 
is also studied, as shown in panel (c) in Fig. ([3]). In- 
creasing the concentration of impurities, we find that the 
peak value decreases. Although the impurities tend to 
suppress the spin Hall conductivity, as would still be de- 
tectable. For a two-dimensional electron gas with its mo- 
bility of order lO^cm^/Vs which is in an experimentally 
realizable regime, the peak value of as is around c/Stt. 
Thus the spin Hall conductivity for a bilayer electron sys- 
tem does not vanish and is expected to be measured in 
samples with high mobility. 

We discuss possible experiments to detect the magni- 
fication of the spin Hall effect. Our proposal is based 
on the fact that a spin-polarized electric current (means 
the existence of spin current) in the presence of the SOC 
can induce different charge populations at the laterals 
and hence a Hall voltage can be detected [l3| . Since the 
induced Hall voltage is in proportional to the spin Hall 
conductivity, its magnitude is greatly enhanced near the 
turning point Pt in the coupled bilayer electron gas. As 
the tunnelling strength can be tuned by the gate volt- 
age, we therefore suggest experimentally detect an enor- 
mously magnified Hall voltage by tuning the tunnelling 
strength to be near f3t in the bilayer system (see Fig.|4]). 
Additionally, the sign changes of spin Hall conductivity 
across f3t also make the coupled bilayer system a candi- 
date for fabricating possible logical gates. 
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FIG. 3: (color online) Spin conductivities in each layer crf(b) 
and the whole system CTs are plotted with parameters r = 
6.6ns, Q+ = 0.55 x 10""eVm, and a_ = 0.45 x 10""eVm in 
panel (a) while q- — 0.45 x 10~^''eVm in panel (b). Clearly, 
sharp cups show up around the turning point /3t. Panel (c) 
is the plot of as with different momentum relaxation times: 
ri = 6.6ns, T2 — 2.1ns, T3 — 0.66ns while the other parameters 
are the same as in panel (b). 



SOC in each layer may be different and the tunnelling 
between the two layers occurs. We gave natural defi- 
nitions of the spin density and the spin current density 
in each layer and derived the corresponding "continuity- 
like" equations. Based on the calculations in Heiscnberg 
representation, we obtained the spin current. The curves 
of the spin Hall conductivities in each layer exhibit sharp 




FIG. 4: (color online) A proposed experimental scheme to 
detect a magnification of the spin Hall effect. By injecting 
spin current (charge current with spin mostly polarized) into 
a bilayer spin Hall bar and connecting the voltmeter to mea- 
sure the transverse voltage, a greatly enhanced Hall voltage 
is expected to be observed near the turning point /3t. 



cusps around the turning point and the peak values have 
signs changed across this point. We also investigated the 
influence of impurities on the spin Hall conductivity. We 
found that an arbitrarily small concentration of nonmag- 
netic impurities do not suppress the spin Hall conductiv- 
ity to zero in a bilayer system, which is quite different 
from the case in the single layer system. The opposite 
signs of the conductivities in the absence of impurities be- 
come the same in presence of impurities. Making use of 
these features, we proposed a possible experiment to de- 
tect a magnified spin Hall effect by direct electronic mea- 
surement. The sign-change property may also be used in 
designing certain logical gates. 
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VI. SUMMARY 

APPENDIX A: EXPRESSIONS FOR SOME 

We investigated the properties of the spin transport COEFFICIENTS AND THE MATRICES 

in a coupled bilayer system where the strength of the 



The coefficient matrices in the linear Eqs. ([8]) are written as 
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The normalization coefficients for the eigenvectors in Eqs. ([6]) read 
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The coefficients in the Eq. ([TOl) are written as 
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The total spin current for each layer in a realistic sample of size x Ly is given by 
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The above results are derived by assuming that the special point fcj is far away from the Fermi momenta. When 
kpz < kt = kpi = kp2 < kpA, the spin current produced by the state with kt is given by ?i^(2a^ — a'^)/4ma^kt in 
unit of eE/An. 

The matrix elements J12, J14, J34 of Jy can be solved from the Ea. ([T6)) . which is in fact a task of solving a set of 
linear equations 
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The solutions are given by 
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where vp is the Fermi velocity. 
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